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Abstract 

We study the density of states method to explore the phase diagram of the 
chiral transition on the temperature and quark chemical potential plane. Four 
quark flavors are used in the analysis. Though the method is quite expensive 
small lattices show an indication for a triple-point connecting three different 
phases on the phase diagram. 

1 Introduction 

To clarify the phase diagram of QCD and thus the nature of matter under extreme conditions is one of the 
most interesting and fundamental tasks of high energy physics. Lattice QCD has been shown to provide 
important and reliable information from first principals on QCD at zero density. However, Lattice QCD 
at finite densities has been harmed by the complex action problem ever since its inception. For /r > 0 the 
determinant of the fermion matrix (detM) becomes complex. Standard Monte Carlo techniques using 
importance sampling are thus no longer applicable when calculating observables in the grand canonical 
ensemble according to the partition function 

ZgM = J DC/detM[U](^)exp{-SG[U]}. (1) 

Recently many different methods have been developed to circumvent the complex action problem for 
small /r/T CHS- For a recent overview see also Ei]. 

2 Formulation of the method 

A very general formulation of the DOS method is the following: One exposed parameter ((j)) is fixed. The 
expectation value of a thermodynamic observable (O), according to the usual grand canonical partition 
function O, can be recovered by the integral 

< O >= y # {Of{U))^p{(^) ^JdcP {f{u))^p{(^) (2) 

where the density of states (p) is given by the constrained partition function: 

p{x) = Z^{x) = J VUg{U)6{4> - x). (3) 

With ( )^ we denote the expectation value with respect to the constrained partition function. In addition, 
the product of the weight functions /, g has to give the correct measure of Zgc- fd = detM exp{—S q}- 
This idea of reordering the partition functions is rather old and was used in many different cases Ellliii 
The advantages of this additional integration becomes clear, when choosing 4> = P and g{U) = 1. In this 
case p{(j)) is independent of all simulation parameters. The observable can be calculated as a function 
of all values of the lattice coupling [3. If one has stored all eigenvalues of the fermion matrix for all 
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configurations, the observable can also be calculated as a function of quark mass (m) and number of 
flavors||5(] (Nf). In this work we chose 

(j) = P and g = |detM| exp{—S'c}, / = exp{i6l}. (4) 


In other words we constrain the plaquette and perform simulations with measure g. In practice, we 
replace the delta function in Equation (H by a sharply peaked potential [ihS. The constrained partition 
function for fixed values of the plaquette expectation value can then be written as 

p{x) ^ JVU g{U) exp {-V{x)} , (5) 

where exp{—1/(x)} is a Gaussian potential with 


y (x) = {x - p)"^. (6) 

We obtain the density of states {p{x)) by the fluctuations of the actual plaquette P around the constraint 
value X. The fluctuation dissipation theorem gives 


u 

— hr p{x) =< X — P >3 
ax 


(V) 


Before performing the integrals in Equation © we compute from an ensemble generated at {pQ,Po): 


{Of{U))Jp„f3) = {Of{U)R{p,po,P,Po))J{R{pf^o,P,Po)),, ( 8 ) 

(/([/)), {p, f3) = {f{U)R{p, /TO, /?, /3o)). / {R{p fio, P, /3o)). , (9) 

-^lnp{x,p,P) = {{x - P)R{p,po,P,Po))^. (10) 

Here R is given by the quotient of the measure g at the point {p, /3) and at the simulation point {pQ, Pq), 

R{p,Po,P,Po) = g{p,P)/g{po,Po) = - Sg{Po)}- (H) 


Having calculated the expressions - dlOb . we are able to extrapolate the expectation value of the observ¬ 
able © to any point {p, P) in a small region around the simulation point {p^, Pq). Eor any evaluation of 
(O) {p, P), we numerically perform the integrals in Equation ©. We also combine the data from several 
simulation points to interpolate between them. 


3 Simulations with constrained Plaquette 


The value we want to constrain is the expectation value of the global plaquette, which is given on every 
gauge configuration by the sum over all lattice points (y) and directions (pu) of the local plaquette (y) 
and its adjoint P},^{y), 


1 r. 


^ = E E « TrP^,(y)+TrP/,(y) 


y i<At<i/<4 ® 

Since the plaquette is also the main part of the gauge action. 


( 12 ) 


Sg 


-/^E E 


TrP^i.(x) + TrPt^(x) 



(13) 


the additional potential V can be easily introduced in the hybrid Monte Carlo update procedure of the 
hybrid-R algorithm f7|. After calculating the equation of motion for the link variables Uf_i{y), we find 
for the gauge part of the force 


iHfiiy) = 


^U,{y)T,{y) 


1 + 


7(x — P) 

P 


(14) 


J TA 


Here the subscript TA indicates the traceless anti-Hermitian part of the matrix. We see that in each 
molecular dynamical step the measurement of the plaquette is required. However, the only modification 
in the gauge force is the factor in round brackets. 
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Fig. 1: (a) Results for Simulations at /3 = 4.98, ^ = 0.3, A = 0.02, n/ = 4, am = 0.05, and number of lattice points: 4“^. 
Shown are the density of states p(x), the phase factor (cos(d)), and their product, (b) Results for Simulations at /3 = 5.1, 
A = 0.01, Uf = 4, am = 0.05, and number of lattice points: 6'^. Shown is the suppression from the complex phase of the 
fermion determinant {cos{6)) for different chemical potentials. 


4 Simulation details and the strength of the sign problem 

Simulations have been performed with staggered fermions and Nj = A. We chose 9 differed points in the 
(/9, pi)-plane for the 4^ lattice and 8 points for the 6^ lattice. On each of these points we did simulations 
with 20-40 constrained plaquette values, all with quark mass am = 0.05. Further simulations has been 
done with (/?, /r) = (5.1, 0.3) on the 6^ x 8 lattice for am = 0.05 and am = 0.03. In order to calculate 
the plaquette expectation value, or its susceptibility, one has to perform the following integrals: 

{P) = J dx xp{x) {cos{9))^ , = J dx x‘^p{x) {cos{9))^ . (15) 

Thus the functions p{x) and {cos{9))^ have to be known quite precisely. We plot both functions in 
Figure^a). The transition is signaled in the double peak structure of p{x). The phase factor {cos{9))^ 
suppresses the peak of p(x) at smaller plaquette values, which results in a shift of the critical temperature 
to smaller values, in comparison with the phase quenched theory. In Figure Q^b) we show the phase 
factor for different chemical potentials. With increasing chemical potential the phase factor becomes 
compatible with zero within errors. In fact, its average value becomes as low as cos{9) ~ 0.005. There 
exist however a small interval around P ~ 2.85, where the phase factor stays finite. In this way, the 
Plaquette expectation values is strongly altered by the phase factor. Figure Q^b) demonstrates also the 
advantage of he DOS method over the other approaches of lattice QCD to finite density. Using the DOS 
method one is able to do simulations at directly those Plaquette values which are relevant at finite density. 
This solves the so called overlap problem of the reweighting approach. Furthermore we have checked in 
f8.l|, that results with have been obtained within the framework of the DOS method agree very well with 
earlier results from the multi-reweighting approach. 

5 The Plaquette expectation value and the phases diagram 

Performing the integration in Eq. O numerically, we calculate the plaquette expectation values as 
shown in Fig.|2l At chemical potentials p<0.36, the plaquette signals the QCD transition through a rapid 
crossover from a low temperature phase of<P>~2.9toa high temperature phase of < P >~ 3.1. 
For ^>0.36 the plaquette expectation value at small temperatures drops to < P >~ 2.85. This new 
low temperature phase of the plaquette at high chemical potentials is caused by the fermion determinant. 
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Fig. 2: Results for Simulations at /3 = 5.1, A = 0.01, Uf = A, am = 0.05, and number of lattice points: 6"*. Shown is: (a) the 
Plaquette expectation value as a function of the coupling /3 for different chemical potentials and (b) the plaquette expectation 
value at fixed coupling, as a function of the chemical potential. 




Fig. 3: The phase diagram in physical (a), and the quark number density at constant temperature T = 143 MeV lattice), 
T = 124 MeV (6"* lattice) and T = 93 MeV (6^ x 8 lattice). 

As on can see in Figure ^b) the region around P ~ 2.85 is the region which is less suppressed by the 
phase factor. Another interesting observation is that the critical coupling, which is decreasing in /r for 
/r < 0.36 starts to increase for /r > 0.36. The plaquette expectation value thus suggests the existence 
of three different phases in the (r,/i)-diagram with a triple point, where all those phases coincide. In 
Figure Ha) we show the phase diagram in physical units. The phase boundaries were determined by 
calculating the peaks in the plaquette susceptibility. Note, that we make no statement about the order of 
the transition lines. To determine the order of the transition one has to perform finite a finite-size-scaling 
analysis. 

The scale was set by the Sommer radius rp, measured on a 10^ x 20 lattice. The triple point is 
located around /i*’’' 300 MeV, however its temperature (T*'"') decreases from ^ 148MeV on the 

4^ lattice to T*''' 137MeV on the 6^ lattice. This shift reflects the relatively large cut-off effects one 
faces, with standard staggered fermions and temporal extents of 4 and 6. 

Also shown in Figure |3l a) are points from simulations with quark mass am = 0.03. The phase 
boundary turned out to be — within our statistical uncertainties — independent of the the mass. 
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6 The quark number density 

To reveal the properties of the new phase located in the lower right corner of the phase diagram, we 
calculated the quark number density, at constant coupling /? and at constant temperature respectively. To 
obtain the density Uq we perform the following integration 


d In detM \ 
d(a/x) / 



/din detM 


p{x) 


(16) 


The thermodynamic quantity Uq are given as usual by 

1 / din detM \ 

a^N^Nt \ d(o/i) / 

In Figure |3b) we show the baryon number density, which is related to the quark number density by 
ns = nq/3. The results are plotted in physical units and correspond to a constant temperature of 
T « 143 MeV (4^ lattice), T « 124 MeV (6^ lattice) and T 93 MeV (6^ x 8 lattice). In order 
to divide out the leading order cut-off effect, we multiply we have multiplied the data with the factor 
c = SB{Nt)/SB, which is the Stefan-Boltzmann value of a free lattice gas of quarks at a given value 
of Nt, divided by its continuum Stefan-Boltzmann value. At the same value of the chemical potential 
where we find also a peak in the susceptibility of the plaquette (pc), we see a sudden rise in the baryon 
number density. Thus for /r > we enter a phase of dense matter. The transition occurs at a density of 
(2 — 3) X UN, where un denotes nuclear matter density. Above the transition, the density reaches values 
of (10 — 20) X UN- Quite similar results have been obtained recently by simulations in the canonical 
ensemble 191. 
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